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rt’e indieste ‘now cq extend th.a SAseh algorithm to handle e class cf 
spatial functions defined in terms of ir.cegrals. 1-iost of the 
integration Ttschinery for this c^ass of funeiiotlS is sitlilar io tha 
rrachir.sfy in the algorithn -whitb handles ibgAPithnSn A progran 
embodying lurch cf the extended integration algbritlha ofiS baap. 
urit les . . It w*s used to check a table of integrals and it snetieded 
in finding sosa niaprints in it. 


The purpose of this ia ta peine cue Cloe c itil Larray Cor cKc 

purposes of integration of the logarithm fiLattica end a CinsS of special 


i unctions defined in Leras of integrals,. Xeae cent 
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that is* the logarithm can be defined as an integral. ftote also that we do 
not res trite ourselves eo a single fraction log (:&)., bnt tet.C fdcr the family 
Or functions logipi^x)). Here we Snail consider the class of special func¬ 
tions F defined 'by 
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where £ is a eons taut t and g(t) can bn obtained by rational operations upon tilt 
Pitiable t by exponential and logarithmic terns involving t. The function 
g may niso tonfain spatial function terms not involving F„ but it mist not 
33 ifitsgrabl* in te-rcs of the allowed functions. This class of special 
functions includes. 
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traditional inccgratlOA methods {a.g., in tog ho cion by pcrta , scba titacios 
of variables) because it is i decision procodL.r l . This means that the 
procedure should be able to ob;cin the integral id in is expressible in. 
closed foro. cr to shew -cEliC the intasra.1 capiat be so caressed. 1’ro.dit,£anal 
^azitods cannot make 5ucU d do-sis ion end are parties; la rly weak if tba 
iat*£re.r,d involves division in a nontrivial way. Or. the other hind cht 
Sis eh algorithm u a shall Oiaploy is weak in that it does not handle, algebraic 
s:-rpz ess ions (e.g,.,, /x). However, such, a restriction has been reeved by 
?.iach in extensions Cb ins algorithm described in [3,4] . 

Oar description cz Ena steps necessary to attend the Risch algorithm 
will c Sneer, trace on the paraicnLar case of the error function. Tfe have 
ex tended the current state of our Implementation of the Risen a Igor i ilia to 
handle the error function. At the present tlhie c«r ictpl*i&Miiition of the 
Ria to d Ijcrithm contains £ hSirly c-Omplcte imp Lensnta Sica of the logarithm 
case and a partial implementation of the exponential case. Details ftbcuS. trie 


the t 


«« taxes £a our implementation are to hi found in [ 2 ]. 

The first thing one hbtitas in attempting an axESnSion is that the 
educEith of Special functions will cams a a modification to be jftade to 
LiQ-aVilie theorem which is basic, to the Risen algorithm, "he TLiouviila 
rtm Says, roughLy speaking, chat the integral of an expression Z it the 
Or Oh SnprashiCn which ls rational in. the logarithmic end exponential 
£ aliitiy in E nd constant multiples of new iogtriihn.it tarns, "he 
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afid logo except constant or Lo^s occurring In a a^r., t 

Ihn generalised versi&n of □'■-2 li&vViHc Esieorfci •jcya ihit coils tint 
ffrpltiplas of special functions can bo ir,eluded in dse integral as wall jo 
esnstent Bsultiplee of logs-. Eti any particular sit-amior; nr 2 ■j&j 16 be loai- 
ftng for integrals invoLving a. given see pi special functions only, say, 
e-zvo-j fiiiictioiLS and exponentitl integral:.. la. even s. case the nraatfa iiatd 
vEfEien e£ ihe Liouvili* cheor&si pt 1 l* u JLiU ikL.1T I_Z 1 Z.S ZO bLirG j_CHTT i. Oil i2~liI: 
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wL'i'sra w y Vj s u i( end t, axe rations,! exomssiers in die exponential* log- 

** 4 11 

eriimic and special inncti&n Cfid=S already in the integrand* and die 
c f , d £ j and are cans He at 1 . 

IP 13 isipQxcaat to note hot-j constant mil. tip Lee or special ft ncGiCnS 
arise in t.t i£teg?£L.+ ■ 2£ia situation in die log case is cults e imp La —'' 
das integrand PUS I hi of the feta 
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The Other a*W log terus {i,£< , not, involving the siiin variable) arise 
through, tha recursive use of the Integra Z L on process in the polynomial 
parti &u*h tirhiS must also be obtained iron rational toms at some 
in the recursion. 

The lem Of the integreod leading to the error function term is 


.1 



ha set ot this form It is clear tear the seta variable as the poles in she 
Integra Eton pros ass when the error function arose was an ■aicpouential. Further 
mote, the term, arose when integrating tee polynomial part of the expression 
or rational zeros of the form 


f x 

<0 


^Tr.e cotter terras are erected lint po^ynem.i&fc terms in the iagO-ritan, whan 
is the mein variable),, 

ha shall now ir.iicrze how the algorithm ter. he absented to handle the 
evrer function. ■„*; u£ suppose hot $-i ,^(;;}) is the iacin variehle in the 
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wir.efa A" find 3 poXjtVOCllali Ln erf (u(kJ> obttinfid ny a ainplft c&spnl 
baaed c“. ills jijcd oi $ find 5 J * For i=lj the denari ft sc or Le eonpletely 
=istored sr.d we proceed fi 5 in ?2j p a f< ft d iosarithaic terms h 
For ass.apla d consider trie integral 
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h jL;i j-E;ccar earn need no; be ticcored since the des.raift&SOS is Utiar: 
:f OiJ * This tarn yields 
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Sow let os consider the p dyne:;, in’. pnrt. 5-oppose t;« iiXpifea&is 
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where the A^ do not involve er£<U>* 2heh the integral, i£ it twites, aye £ b& 
Or the fora 


■ n+: er/ l+l (u) * i erf 3 (-4} -)-+ + +-r , 


vjL’-ere ?'■’.? 3. do not involve erx(u) , aad only 3„ a:sy cor.ceLa. new lo-saxitlrsie 
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of crrof-fuafitlftn teras, 

Wa cun ohinin values for the 1. (nr decide the nonlategrabi11cy of the 

L 

iiioftsaisn in Logori thaic, exponentiali of error “function terma) by a process 
similar to the One described in [2] , An erupie which uses this Method is 
given below. Consider 
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o- a + V = 1 


By in tests, ting (trivially on the loft-Sind, by a recursive, call EO Ena 
in zc.grs ti On pric&dsro Or. Lr.^ ri^bt), kb Obtain 
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integration yields 
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(Saa below for a discussion of such th integration stop). 
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Subs c fitting th-aSe results into the genere1 foti of the integral, we 
obtain 


x art* 1 fa) -r —^ e X erf fa) - / ™ ertfa /l) + b, 


It thus Altars that when the main variable is erf fa fa)} the algorithm 
la ejite a ini her to the _og ease as ve have claimed. iTna algorithm must also 
he modified in the exponential ea&i, TCa have to permit eons tent uultipLeg 
Oi error-fanetiOn tarns to he fan faded fa the.solution. This follows from 
cor previous discussion, Ihit fa* if the main variable ia a~ v "'| then 
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Hence 


A h B c + klsi r -p p 1 


£.rs given, A „ Up c*.c is. Our fra j-act is *.a find L and z , Vfi&a these 
, , 2 

v.Tcra -p " ku. 3 is rtstriecai to beir .3 rational ia the logarithmic p 
Enpecen..J-Sij and o^iOr - rune,t ion tama which A contn ilia . ISfrtefrvar t e in-ce ua 
do not ailcw algebraic empress iocs &Kii*S for Elgebraie tfrr .5 rants, ut chall 
restrict p accordingly, Rota that -Isu. cin. be skprsEsed es £ rs-Sionai. azpras- 
Sivn in the variable o£ integration and any additional defined Variables 
raprasantfag io$5 , fcfipOr.anitials and error-functions - He want to know if -ku 
ls a peri;-;tt Squire^ This can be dntided by factoring. If ~kc is a perfect 
sc tiara | let p be a sc care “&QS of it. Otherwise ignore the error fjneciQ^ 
terns and proceed AS uSudl„ 

Re orrinsing the equation above, we obtain 


R 1 -r hoc 1 = A - 
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7n 


. 'a can new ns a tee -itehinary of the Risen a Lgorithst to solus for 3 end c, 
-- may bistst. A sotorien for E and o coat be unicue for if ic ware not it 
tOu-c rnply tent an error-function rorvi is □ Ljpabraitally dependant on an 
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7cr £iOi special function cat mm t perform an analysis of the ways in 
which corns cans t,u 11 ip loi df an instance of it ten catat inco the. integral. 
T/.is involves some pattern recognition, oat la soema difficult to cbiSin a 
a imp la algorithm to perfopa die £c£k in genera!* However, for itch special 
function, definable by an Integra L that we h^vo encountered, the analysis is 
quite s imp 1.0 [see 5] , 


We must also know in whet ways the ipooial 
previous functions, The algorithm assumes the 
the terms. Hence t we must not let a Lge'cra ica 11 
appear in fits inte^ranci, An emampie of seen a 
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{ log (h 


function interacts with 

al^ebrtit independence of al 
"j dependant omprescions 
dependence is gives below- 



'ihere may also be algebraic relationships between different instance® of 
z r \m strae function log{e.b} * log e + log b + Zksri}* in these cases 

Wi presently have to Tely on previous analyses of the simpLifitstlon rules 
Sts otic ted with the special functions tnder consideration. Divin^g the 
simplification rules applicable to soy gives special function reiativa 
to a given set of fane cions sasc.s vary difficult it present. 

It should also ba clear that tbs fact that an expression is nor integrable 
it terras, cf existing rur.crier.s doc; sot necessarily miht its Ir.iegrel e 
Suitail* othdidfeM for a Special function, 70S 4aamplft,' net* rhec 
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because of tills re lutiouship the integral W. Elia Loft is z bad candidate 
ior £ special function. A lit fit GOr.b id Station will reveal several restric¬ 
tions flu what cna (Sight Gill "primitive" integral ■ fogtsicclft!} ipGcifcl 
funsttonfl to these defined by such pi'idiciVi La Cssft La fill -.take the eaten- 
sion e good deal c-aiicr* 

In order to dtfcug tha pxogrc.7.1,, we uaed our version of the Sisth al^orithin, 

er Loaded to hir.d ia the error iuacuiOb to thick A line table of Integra is Or 
the «sor function compiled by \-j. 2. Knurcr ia 19>3 1$]* (A mere recent 
tGbIe a emphasising definite integrals is {7 ]), The program unfortunately 
cannot check the correctness Of the entire set of integrals because acc?.e 
entries involve algebraic expressions such as /* and others involve dafinite 
ntegrals + Ina program has so fax yr, covered a sispri^e ih 


t 2 . fi< 2 l\ 2 ^ x 
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2 2 
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chore the under lined express loft was 2 /"li. 

.-. potentially controversial disagreement b£Ci*ee3 the program and the 
table occurred in 


! :-i^ erf (ax -T- b) dx * -V sri (ax +"3} s^ K (2 - ipn + p^x^) 

J -I 5 

2 

^-r 9 

_ rax - nh + 2a - £e . - (ex^b) - px 

o 2 e 

*p A 




GiIT - \^m%. T S' 



a 








«._QV%jVi31“ , _TC L" ‘ i _.i^ B ... n . 


- i ■ 


no ii£- ‘i <L 


!.« &Z.L-Z. wJal t'-'-.i 

iif.dZuC LQii VIW - jJ C. O.J.d liVv t ' - U vH v vfilv --■_ Z." 

-'■'-' - t '■-' Ciltiii- 1 - -■. - '■-■ IFjh^Iluv iHu'.Zi ■ V.V-- 

lO enter LM0 arena uith a coiTti.ij.iZiitior: or (ii.ii L Civ oz the proroscd jc 1-j f 

at wisSs CP other solution of zhe.tr cihooa i&g. 

ftauraf recently ir,£asrqei tai that he tpanz at least zorcy hours ciucU&Rg 
by h&ud the solution to the ab ova ’ircblezi before he gave up Of, if. Inc ?DI?-5- 
i.; i'll spent eight nicL-tcs on the prop Ito; i IJi btliove mat the nucnir.e's tin. 
cur. be reduced to n raw seconds by rating changes to the eimpli£ier as ad Lti 
shs program- Currently the a imp II tier re rep recer.cs expressions quite iracuer 
due to a real" that ftCw logari chitio or error function tEir.is were obtained 
in the dower levels *£ tire recursion. The algebraic manipulations in the 
program ers ttostiy performed using a rational function rentesentacion, but 
this rentesentetion does not elbow one to introduce new variables, ^ner. a 
new variable appears , the old expressions xist he rerepresented to account 
£ct its exisioz.ee, ForcunataLy the sinplifiar's scats are usually ground- 
Lena, ana appropriate medizioati bus to it should yield considerable savings 
in ei:e cution t ime . 

Sisth has pointed out to us shss 0® trowski [&J gives aa algori;h?i which 
is* ir. the case where no algebraic expressions appear, an extension a; the 
His eh aigorirlM for the leg cate- Cstrewski deals with integrals of rational 
entreesions involving a function defined in terns of an integral, [iowever, 
he considers only a single function (e.g. * log(x) or eif(]{ ™ 1) ) rather 
taan a zaz.tly oz ZuZicZzons■ nenae za caars wttn. the case when a -.eg or 
soatza* zuucczon ts cue main vartaale oz zr.tagractot j but coat not nnndIs 
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.-,a ft have seenj the problems that rautia tr.seIvan an such cases are ssstiU- 
tially 3±K‘.pti£ica,,v.i .H. 
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